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Abstract 

In this note we demonstrate that the algebra associated with coordi- 
nate transformations might contain the origins of a scalar field that can 
behave as an inflaton and/or a source for dark energy. We will call this 
particular scalar field the diffeomorphism scalar field. In one dimension, 
the algebra of coordinate transformations is the Virasoro algebra while the 
algebra of gauge transformations is the Kac-Moody algebra. An interest- 
ing representation of these algebras corresponds to certain field theories 
that have meaning in any dimension. In particular the so called Kac- 
Moody sector corresponds to Yang-Mills theories and the Virasoro sector 
corresponds to the diffeomorphism field theory that contains the scalar 
field and a rank-two symmetric, traceless tensor. We will focus on the 
contributions of the diffeomorphism scalar field to cosmology. We show 
that this scalar field can, qualitatively, act as a phantom dark energy, an 
inflaton, a dark matter source, and the cosmological constant A. 
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1 Introduction 



Recently experimental cosmology has brought new challenges to the theoretical 
understanding of the universe. Until now, the standard hot big-bang cosmology 
has been very successful in explaining four important observed facts: (1) the 
expansion of Universe, (2) the origin of the cosmic background radiation, (3) 
the nucleosynthesis of the light elements, (4) the formation of galaxies and the 
large scale structure of Universe. However, recent discoveries from Supernova 
Cosmology Project and WMAP strongly suggest that the standard 
cosmology needs to be revised. This new standard cosmology must now explain 

o Flat and accelerating Universe 

o Rapid expansion of Universe in early period (Inflation) 

o Composition of Universe: 73% dark energy; 23% dark matter; 4% baryons 

The theoretical search for this new standard cosmology is very active and there 
is an overwhelming amount of literature published where most of them have the 
same theme: scalar fields. In this work we ask if there is a guiding principle that 
might account for the presence of new fields in cosmology. We posit that there is 
a symmetry that can accompany general relativity in providing an explanation 
for some of the new features found in cosmology. This symmetry is not new and 
is already exploited in general relativity through general coordinate covariance. 
However we will realize the symmetry using an approach used in string theory. 
Since string theories are two dimensional field theories, direct contributions to 
gravity coming from curvature and the Einstein- Hilbert action are diminished 
because of the low dimension. In fact in two dimensions the Einstein-Hilbert 
action can at most distinguish different topological structures. Our starting 
principle will be the algebra of Lie derivatives which exists in any dimension 
including one dimension. We will show that this algebra admits a field theoretic 
representation analogous to the Yang-Mills gauge potentials that arise from the 
algebra of gauge transformations. 

In one dimension, this field theory comes from the coadjoint representation 
|15ll6l35fl^ of the one dimensional algebra of Lie derivatives called the Virasoro 
algebra. In string theory this representation is used to understand the geometric 
origin of gravitational anomalies. As we will review here, one can understand 
the origins of two dimensional gauge and gravitational anomalies through the 
Geometric Action |4l31l20l21l26j as well as establishes the Lagrangian for the 
fields living in the coadjoint representation called the Transverse Action |6I8| . 
For this work, the importance feature of the Transverse Action as compared to 
the Geometric Action is that the Transverse Action can exist in any dimension. 
We will review the construction of the Transverse Action with both the algebra 
of gauge and coordinate transformations. We show that there is a naturally 
occurring scalar fields called the diffeomorphism scalar field (diff scalar) that 
arises naturally from the algebra of associated with coordinate transformations. 

As a matter of completeness we will keep this paper somewhat self contained 
and outline of this paper as follows. We will first review the Robertson- Walker 
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cosmology and as well as four of the more popular scalar field theories discussed 
in the literature. This will introduce our notation to the reader and put this 
work in a more familiar setting. In the next section we discuss the method we 
use at arriving at the Lagrangian for the diff scalar. We are then able to discuss 
the cosmological implications of the diffeomorphism scalar field by combining 
the action for the diff scalar with the Einstein-Hilbert action in the next section. 
We show that the diff scalar field can, qualitatively, act as phantom dark energy, 
an inflaton, a dark matter source, and the cosmological constant A. 



2 Standard Cosmology Review 

The standard cosmology is based on two postulates that Universe is homoge- 
neous and isotropic at large scales |24I33I12| . These postulates are strongly 
supported by the cosmic microwave background (CMB) radiation coming from 
various regions of sky: the temperature of the CMB radiation from different 
parts of sky are almost identical. These postulates are called the cosmological 
principle. With proper choice of coordinates these symmetry postulates thus 
allow us to write the metric in the form of Robertson- Walker metric: 

ds^ = -dt^+ (t) I + {de'' + sin^ 6 dc^^) | (1) 

Here we used the convention: Ti — c — 1. The function a (t) is called the scale 
factor of the universe {cosmic scale factor). Rescaling this factor appropriately, 
the constant k can be always taken to be 0, +f , or —1: 



k = +1 : closed universe with positive spatial curvature 

fc = : flat universe with zero spatial curvature 

fc = — 1 : open universe with negative spatial curvature 

Since current observations seems to prefers a flat universe, we will adopt 
k = throughout this paper. The cosmological principle also demands that all 
cosmic tensors are maximally form-invariant with respect to the spatial coor- 
dinates. As a consequence of this constraint the most general form-invariant 
tensor for the energy-momentum tensor is given by: 

Ton^pit), Toi=0, T,j=p{t)g,j, z,j = 1,2,3 (2) 

where p (t) and p (t) are arbitrary functions that depend only on time. Rewriting 
the energy-momentum tensor as 

T^_,„ = ip+p)U^U^+ pg^^ (3) 

where is a "four- velocity vector" 

C/° = 1, W = {i^ 1,2,3) (4) 
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we identify p (t) and p (t) as with energy density and pressure of a perfect fluid. 

Using this form of the energy-momentum tensor and the Robertson- Walker 
metric, we obtain the Friedmann equations from Einstein's field equations: 



^ - 3M2 +3 aHt) ^ ^ 

where the Hubble parameter H{t) is defined as 

and 

Mpi = reduced Planck mass = (SttG)"^/^ . (8) 
From the energy-momentum conservation condition V^T^^ = we find 

p + 3H{p + p)=0. (9) 

This equation and the Friedmann equations are not independent from each 
other. The equation of state, the pressure-to-energy density ratio, is defined as 

w=-. (10) 
P 

Typical values of w are 

Radiation: w = ^ (11) 

Matter: w = (12) 

Cosmological constant: w = —1. (13) 

The condition for the accelerating expansion of Universe can be given in 
several ways: 

the cosmic scale factor is accelerating : a > 

d fH^^{t)\ 

the comoving Hubble length is decreasing : — I — — ) < (14) 

dt \ a[t) J 

negative Pressure : p{t) + 3p {t) < 0. 
The last condition can be also expressed by using the equation of state: 

w<-l (15) 
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2.1 Single-Field Inflation Model 



The inflation models assume an era of inflation during the Big Bang, in which 
the energy density of the Universe was dominated by the potential energy of the 
scalar fields (inflatons). The inflation is defined as any epoch during which the 
conditions Eq^l are satisfied and used for the events in early Universe. The 
negative pressure condition for the inflation model motivates one to introduce 
a scalar field (inflaton) with its Lagrangian density given by 

^0 = ^5^"^ (9,0) - y (0) . (16) 

The term V {(p) is the potential of the scalar field and the choice of its form 
characterizes the types of inflationary models. Using the above Lagrangian 
density we have 

T^- = (9''0) (d'^^b) - g^" (9V) (5a 0) + V (</.)] (17) 

or 

p = + 1 (v,^0)2 + V {^) (18) 

and 

P=h^'-l (V.0)' ~V{^). (19) 
The equation of motion for the scalar held is given by 

^ + 3i70--^V^0 = -|^. (20) 



2.2 Quintessence 

In quintessence Q is a scalar field that is introduced to explain the acceler- 
ating Universe. Quintessence is a dynamical field and generally has a time- 
varying, spatially inhomogeneous negative pressure [H]. The Lagrangian for the 
quintessence is given by 

£ = ^d^Qd^'Q ~ViQ). (21) 

The energy density and pressure of quintessence are, then, respectively, 

p^lQ^ + V{Q), P=\Q^-V{Q). (22) 

It is easy to see that the ratio of the pressure to the energy density, w, is —1 < 
w <Q for quintessence. A particular class of quintessence models have a tracker 
behavior. In these models the density of quintessence closely tracks the radiation 
density until the equality between matter and radiation is established. Once the 
matter-radiation equality is established the quintessence start to behave as dark 
energy, k-essence and phantom energy are the special cases of quintessence. 
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2.3 k-essence 



The idea of k-essence (or k-inflation) was originally introduced as a possible 
model for inflation [p. In this model an inflationary evolution of the early 
universe is driven by non-quadratic kinetic energy terms starting from rather 
arbitrary initial conditions instead of the potential energy term V {(p). When the 
idea of k-essence is applied to the coincidence problem-'^ , it was shown that 12 
131 k-essence can act as a dynamical attractor at the onset of matter domination 
period and introduce the cosmic acceleration at present time. A purely kinetic 
k-essence model was investigated in |2H1- In this investigation the author showed 
that k-essence in this model can evolve like the sum of a dark matter component 
and a dark energy component if the Lagrangian for the k-essence has a local 
extremum. 

The Lagrangian for a typical k-essence model is given as the pressure of the 
k-essence: 

C^P^V{<t>)F{X), (23) 
where (p is the scalar field representing k-essence, and X is defined as 

X = hd^,4>d^'(j). (24) 

The pressure in this model is given by the Lagrangian itself, while the energy 
density is given by 

p = 1/ (0) [2XFx - F] (25) 

where 

= f . (26) 

Using these pressure and energy density, we have the ratio of the pressure to 
the energy density, w, as 

The behavior of this ratio is determined by the form of V {(f)) and F {X). 



2.4 Phantom Dark Energy 

The phantom dark matter is a special case of quintessence characterized by the 
ratio of the pressure to the energy density, w < —1. Though the dark energy 
with w < —1 will violate all energy conditions [221, the dark energy models 
with w < —1 that are consistent with observed data can be constructed |10j . 
A phenomenological model of phantom dark energy can be constructed to give 
new insight to the coincident problem through the interaction of the phantom 
dark energy with dark energy and dark matter lllj . 



This is the problem related to the precision required for the initial value of the total 
energy density of the Universe so that the spatial curvature of the Universe is almost zero: 
the flatness problem. 
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3 The Origin of the DifFeomorphism Field 



In this work we posit that there exist a scalar field, the difFeomorphism scalar 
field, that has its origins in one of the most primitive symmetries used in physics; 
general coordinate covariance. When coupled to a Robertson- Walker metric, 
this scalar field qualitatively exhibits the behavior of an inflaton, phantom dark 
energy and the cosmological constant |27j . Here we review the construction 
of the Lagrangian for the diffeomorphism field. We will marry the general 
coordinate transformations with gauge transformations to show the analogy 
with the vector potential in Yang-Mills theories. 

Recall that under a general coordinate transformation, tensors transform as 



^"(-') = (^)^''(-)- (29) 

Their infinitesimal counterpart, i.e. a;'" = x"" + where ^ is vector field, 
determine the Lie derivative; 

C^Aaix) = -C{x)''dbAaix) - Ab{x)dai\x) (30) 

C^A'^ix) = -^{x)''dbA^{x)+A\x)dbC{x). (31) 



It is important to note that for tensors the Lie derivative is independent of the 
connection, Vq. Thus the field theory associated with diffeomorphisms will be 
independent of Riemannian curvature. 

The analogue of the above for finite and infinitesimal transformations for 
gauge transformations on a (non-Abelian) vector potential are 

V^ix) = U{x)Va{x)U{x)-^ + iU{x)daU{x)~^ (32) 

and 

5Mx) = I [k{x), Va{x)] - daK{x), (33) 
where U{x) — expiA(a;) is the finite gauge transformation generated by A(x). 

3.1 The Adjoint Representation 

The vector fields (^'') and gauge parameters (A) act on themselves and form 
the adjoint representation. The gauge and coordinate transformations in one 
dimension are particularly interesting because the dual representation of the 
adjoint representation, called the coadjoint representation, reveals the gauge 
and gravitational anomalies that appear in two dimensional field theories and 
string theories. Interestingly enough, the coadjoint representation appears to 
come from a field theory that need not be constrained to two dimensions. In 
the case of the gauge transformations, that field theory is Yang-Mills. 

To see how these field theories comes about, lets begin by marrying the 
one dimensional infinitesimal coordinate transformations (Lie derivatives) and 
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infinitesimal gauge transformations together. This marriage is the semi-direct 
product of the Virasoro algebra [30] , which corresponds to Lie derivatives of one 
dimensional vector fields, and an SU(N) Kac-Moody algebra, which corresponds 
to SU(N) gauge transformations on a circle or a line. This algebra admits central 
extensions (or phases at the group level) so that the basic elements of this algebra 
contains the Lie derivative generated by the one dimensional vector field 
a gauge parameter A(0), and a constant, say a so we can have a three-tuple 
(^;A;a). Then the commutation relations, [[*,*]] for the adjoint representation 
is 

m- A; a), (r;; S; h)]] = (C o ,7; -^'SfcS + ifd,K + [A, S]; {A, S, C, r,}), (34) 
where the central extensions {A, S, ^, 77} is given by 

{A,S,C,r?} = ^ J{AdaE-EdaA)de'' + ^ J(C^a^b'7cV'')de' (35) 

The parameter h,k, and c are arbitrary constants that define the central exten- 
sion. The new vector field (^ o rj)"' is defined through the Lie derivative, 

(e o rjr = C^V^ = -^'dbif + v'db^- (36) 



3.2 The Coadjoint Representation 

Associated with the adjoint representation, we can define its dual by writing 
an invariant product |15I35) . Let (77; S; b) be an element of the adjoint and 
B = {D,A,iJ,) be an element of the algebra's dual. Then the one dimensional 
line integral gives an invariant pairing: 

{{D, A, fi) I {rj; S; b)) = J {Dab{0)C{d) + Ahie)A{0)) dO" + b ti. (37) 

Demanding that the pairing to be invariant under gauge and coordinate 
transformations is equivalent to acting on the pairing with an adjoint element 
and getting zero. Consider the action of the adjoint clement = (^ (6*) , A (9) , a) 
on the pairing. Then the invariance of the pairing 

S^{iD,A,^^)\i7^;E;b))=0 (38) 

and Leibnitz rule implies that 

S^B = {5^D {6) , 6:fA {9) , 0) (39) 

where 

shift in tlic diff coadjoint clement 
. " . 

(5^D {6) = 2e'D + d'c + ^e'" + - Tr^AvV) (40) 

coordinate trans 9^'^9'^ trans 
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and 

shift in the gauge coadjoint clement 



= ££+£A- [AA- AA] + fc^A' . (41) 

coord trans gauge trans 

We have suppressed the tensor indices for the fields Dab and Ab and ' denotes dg. 



3.3 The Geometric Action 

The fields A in the coadjoint representation correspond to the spatial compo- 
nents of a gauge field while D is the space-space component of a rank two tensor 
called the diffeomorphism field. The "/i" is the central extension. Now there 
are two distinct Lagrangians that one can construct using the structure of the 
coadjoint representation. The first is called the Geometric Lagrangian. 

In the geometric Lagrangian or geometric action ^| [23 the coadjoint ele- 
ments serve as background fields while the coordinate and group transformations 
become the dynamical fields. There one considers a two parameter family of 
coordinate and gauge transformations that can act on a fixed element of the 
coadjoint representation, say B = {D [6] ,A{d),fi). The two parameters, say 
A and r, correspond to a natural coordinatization of the orbit of B. The orbit 
of B is defined to simply be all gauge fields and diffeomorphism fields that are 
related to B via gauge and coordinate transformation. 

On the orbit of any coadjoint element, there exist a natural symplectic struc- 
ture corresponding to a rank two antisymmetric tensor |17I18) . When in- 
tegrated on a two-manifold this will give the geometric action. To find the 
geometric action for our case, namely the semi-direct product of the Virasoro 
and Kac-Moody groups, we define two adjoint vectors corresponding to changes 
in directions corresponding to coordinates A and r on the orbit. These adjoint 
vectors are u\ and u-r- Now let g represent the group action of both an arbitrary 
gauge transformation, U{d,X,T), and coordinate transformation, s{9,X,t). We 
denote by Ad{g) the action of these gauge and coordinate transformations of the 
adjoint representation and Ad*{g) the group action on coadjoint elements. Then 
using two adjoint elements corresponding to changes in the A and t directions 
and a coadjoint representation Wg defined by 

= {^x (9), Ax (9), ax) 
Ur = (^r W,A,(^),a,) (42) 
Wg = Ad*{g)w^Ad*{g){D{9),A{9),fi). 

We then pair these elements via {wg \ [ux,Ut]) to define a natural symplectic 
two form in the coordinates of A and r. This is given by 

{wg \ [ux,Ur]) = {Ad* {g)w \ [ux,Ur]) ^ {w \ [Ad [g^^] ux, Ad {g^'^) Ur]) 

(43) 

Concretely we have 

^x{9) = {dxsi9))os-' and Ax (9) = (OxU) o U-\ (44) 
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and we find 

Ad{g~^)^^{e) = s-^o{d,s{0)) = d,e^^ (45) 
Ad{g~^)Ax{e) = U-^o{dxU)^U-^dxU. (46) 

Then in Ea. (|43|l we can write 

[Ad{g~^)ux,Ad{g~^)ur\ = 

[(iti, U-^dxU, ax) , (fe, VdrU, ar)] = (47) 
([i^'l^] . [li^^U-'drU] + [U-'dxU,%^] + [U-'dxU,U-'drU] ,cv + ck) 

where cy and ck are the central terms due to the commutation relations for 
the Virasoro and Kac-Moody algebras, respectively. Choosing c = —h , the 
central term for the Virasoro commutator is given by 

p27T 

cv = ^ dau'v" (48) 
where u = and v = whereas the Kac-Moody commutator gives 

Cfe = ^ / daTr [uv') , (49) 







where u = U ^d\U and v = U ^drll. Hence, after some calculations, the 
geometric action 



Sv($KM= I d\dT{w\ [Ad{g ^)ux,Ad{g ^) w^] ) (50) 
becomes I4llj<jl:^()lljll 



S j rf'eD(0)(§f) 



Metric Coupling 

2¥ 



hjd^e TrA [6) (l^ae {U~^drU)) (51) 
j d30TrA(0)(|^ae {U-^dxU) + [U-^dxU,U-^drU]) 



+ d^9TT{[U-^deU,U-^dxU]U-^drU), 
where the measures d^6 — dOdr and d^6 — dddrdX. 
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The fourth summand is recognized as the Liouville-Polyakov action [22] while 
the fifth and sixth terms are the Wess-Zumino-Novikov-Witten actions [21 ■ In 
the summand called the "Metric Coupling", one sees that the diffeomorphism 
field interacts with the metric through 

where hat is the two dimensional induced metric. This interaction is the first 
ingredient we will need to discuss cosmology. When the trace of Dab becomes 
constant, this term in the Lagrangian will appear as a cosmological constant. 

3.4 Transverse Actions 

Every point on the coadjoint orbit of S = {D, A, ii) corresponds to a group 
element, s{6) and U{9). There are, however, group elements that do not appear 
on the orbit because they do not change B = {D,A,^). These group elements 
are generated by the isotropy algebra of the coadjoint element, B = (D, A, /i). 
An adjoint element T belongs to the isotropy algebra B — {D. A, fi) when 

SjrB (9) = 2e'D + D'e + + - Tr (AA') = (52) 

and 

SyrA{e)= A',e + ^'A- [AA- AA] + fc//A' = 0. (53) 

These adjoint elements are candidates for the conjugate momenta of the fields 
D and A that appear in B = {D,A,^). To see this recall that in Yang-Mills 
theories the finite gauge transformations are given by 

A,ix) ^ Uix)A{x)U-'{x)--{d,Uix))U-'{x) (54) 
E'{x) U{x)E'{x)U^^{x). (55) 

In temporal gauge 

Ao = 0, (56) 

the residual gauge transformation are time-independent gauge transformations. 
The infinitesimal residual gauge transformations are: 

6Aadx) = [A{x),A{x)l^-d,Aa{x) (57) 
SElix) = [Aix),Eix)]:. (58) 

The isotropy algebras are given by setting these equations to zero: 

5aA„(x) - [Aix),A{x)l^-d,Aa{x)^0 (59) 
S^Elix) = [Aix),Eix)l = 0. (60) 

In Ea. (|54|l and Ea. (|55fl . the transformation laws reveal that in one space 
and one time dimensions, the gauge potential is in the coadjoint representation 
while its conjugate momentum, the electric field, is in the adjoint representation. 
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Because Ei is in the adjoint algebra it can gauge transform Ai. The Gauss law 
constraints 

V^El {x) = d^El (x) + [A, (x) , E^ {x)] ^ = 0, (61) 

forbids this from happening. From Eq. H59I) we see that the conjugate momentum 
to the coadjoint element is an element of the isotropy algebra. With this one 
now has a way of finding a Lagrangian that will give dynamics to the fields in B. 
We construct a Lagrangian where one of the field equations reduces to a Gauss 
Law constraint which corresponds to the isotropy equation in two dimensions. 
The Lagrangian can exist in any dimension. For the Kac-Moody sector the 
Lagrangian is the familiar: 

Sym = ^^ld"x Tr (F^, {x) F^^'' (x) ) . (62) 

This approach was applied to the isotropy equation obtained from the Vira- 
soro algebra, 

d^D = 2D'i + D^' + q^'" + /3C' = (63) 

where 

q=^ and /3=^. (64) 

Following the Yang-Mills case, we assume that this equation is obtained after 
enforcing the temporal gauge and on finds the transverse action, 

^ - J d"x^ a (X'™'' D\X^ia + 2X'™'D,,X\„,) (65) 

where X™"'' = -^rj-frnn^ The term proportional to q arises from the central 
extension in the algebra. Algebraically the central extension exists only in one 
dimension. However, in the construction of the transverse action, this contribu- 
tion is seen as coming from a derivative interaction of the connection with the 
diffeomorphism field. We expect this term to be relevant in the early universe 
when there is only time dependence in the system. As we will see later, after 
time evolution, the q dependence of the system becomes negligible. This could 
be the signature that the symmetries in the cosmological principle are broken 
and that spatial dependence becomes relevant. 

A pictorial interpretation of the geometric and transverse action can be 
seen in Fig.l^. The geometric actions is the physics of the coordinates and 
gauge fields (collective coordinates) and is represented by horizontal slices on 
the foliation of the dual of the algebra. These horizontal slices are equivalence 
classes of fields and any physics there would constitute anomalies in the gauge 
and diffeomorphism symmetries. One the other hand, the transverse actions 
move the fields from one equivalence class to another and do not contain spurious 
gauge and coordinate degrees of freedom because of the Gauss law constraints. 
It is for this reason that we say that the transverse actions are transverse to 
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The Relationship Between Geometric Actions 
the Transverse Actions 




T 

Figure 1: The Foliation of the Dual of the Algebra 
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the geometric actions. In n dimensions the mass dimensions for the couphng 
constants and the difF field are: 

[a] = [13] = 0, [q] = -2, and [D^,] = (66) 

To see that the diffeomorphism action recovers Eg . (163(1 in two dimensions, 
we vary the action with respect to the space-time component D^q. Then by 
setting Dio = 0, the field equation becomes 

XlmOg^j^lrn _ d^{X^"^^Oj^u) - a,(X"'OD„„) - q ^A^„^X'"'° = 0. (67) 

In 1 + 1 this corresponds to the isotropy equation found on the coadjoint orbit 
where D corresponds to to space-space component, Dn. The field equations in 
l+l dimensions reduce to 

XD' + 2X'D + q X'" = ^D' + 2$'D + q = 0, (68) 

where the adjoint element ^ corresponds to the conjugate momentum, X = 
X™ ^ of D = Dii. This action is exactly analogous to Yang-Mills theory. 
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4 Cosmology and the DifFeomorphism Field 



We have seen that there are two types of physical action naturally arise from 
the coadjoint representations of the Virasoro algebras (diffeomorphisms) and 
Kac- Moody algebra (affine Lie algebra). One action is the geometric action 
that comes from the coadjoint orbits of algebras. The other action is from the 
phase space that is transverse to the coadjoint orbits: the transverse action. 
The transverse action of the Kac-Moody algebra is the Yang-Mills action in 
two dimensions and it is also valid in higher dimensions. The transverse action 
for the Virasoro algebra can be constructed by following the analogy of the 
construction for the Yang-Mills theory. The action associated with the alge- 
bra of diffeomorphisms is also valid in higher dimensions as its origins are Lie 
Derivatives. 

The Virasoro algebra is the symmetry under coordinate transformation on 
the unit circle, i.e., diffeomorphisms. There an important lesson from string 
theory is that diffeomorphism should play an active role in the theory of grav- 
itation. Indeed, we have already seen important implication of diffeomorphism 
in the geometric action for the Virasoro algebra Eq. [^l] : 

^ I drdeD{9){l^). (69) 

If we rewrite this term in coordinate invariant notion, we have 

J VhDabh''^dx (70) 

where /i"'' is the inverse metric on the two manifold. Our interpretation of the 
diffeomorphism field D is that it is a non-Riemannian contribution to gravita- 
tion: it appears already in two dimensions where the Einstein-Hilbert action 
only offers topological information. It is easy to observe that when the trace 
of Dab takes a constant value, this term appears as a natural source for the 
cosmological constant. Such a term in the action does not depend on the two 
dimensional structure. 

In this section we will apply the transverse action of the Virasoro algebra 
to cosmology. More specifically, we will study the coupling of the trace part of 
the diffeomorphism field D to the Robertson- Walker metric for the flat Universe 
(k — 1) and study the cosmological implications of the diff field. Many of the 
calculations in this section are done by Mathematica with MathTensor package 
j23l37| . The action used to study the cosmological effect of the diff field D has 
two parts: 

5 = 5hb + 5d.« (71) 

where She is the Hilbert-Einstein action without cosmological constant and Sdnf 
is the transverse action for diff field: 



5he = - -r^ / dx* V^)R (x) , x = (f , t) (72) 



15 



and 

- j (qX'^\ViV^X'^a + ^X^a-Xkga) (73) 

We note that a term is from Gauss' law constraint, q and [3 terms are from 
the central extension of the algebra. A term is analogous to Eq[7n| i.e., the 
cosmological term. Since we will set the ordinary cosmological constant A to 
zero, there will be no ambiguity. Again the mass dimensions in n dimensions 
are 

[a] = ^, [/3]=0, M = -2, [X]^^-^, [^M^]-^. (74) 

The equations of motion and the energy momentum tensor for the diff field 
D fj^u are obtained by varying the diff action with respect to D^i, and 
respectively: 

equations of motion: — (75) 

and 

energy-momentum tensor: (JS'diff ^ ^ J d'^^V 9 (^) '^^'^ (^) ^9fj.v- (76) 

The trace of the diff field D^i, is defined through the decomposition 

D^,{x) = :^^Pix)g^, + W^,ix) (77) 

where i/j (x) is a scalar field, W^,^ (x) is a traceless symmetric tensor field, and 
n is the dimension of the spacetime. Since we are interested only in the trace 
part of the diff field 1?^,^, we set the traceless part to zero: W^i, (x) = 0. The 
equations of motion for Z?^,y (x) , reduce to j38| 

A-^-2a(^ + ^)^ix)]v^V,i,{x) (78) 



n n ^ 

-a (\ + - ) V^V (x) V.V (x) - -gVV'^V.V^V' (x) +nX = 0. 
yn"^ n J n 

where is the covariant derivative. The energy-momentum tensor for D^i, is 
reduced to |2H1 



+ 4 



VV^i) {x) VV^P (x) + V^ip {x) V^V^V^V (a;) 
-V^?A [x) V^VV^V (a;) - V'^V [x) V^V^V^^A [x) 
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(79) 



Applying the cosmological principle to the diffeomorphism field, we write, 
in four dimensions, 

it) = y{t)9^.u (80) 

where f{t) = ip {x, t) is a scalar function that only depends on time. Using this 
form of diff field, we find the equation of motion for the diff field as 

(l - I - ^af (t)) f" (t) + [(l - I - ^af it)) r [t) + {t) 

-UH^nt) + UH'f'it) - hrm'H = o 

(81) 

Here we used ' ' ' as a time derivative and H = H (t). The energy density and 
pressure of the diff field are also found, respectively, as 

p,,,(i) = ^(l-^-^af{t))f'{tf + Xf{t) (82) 

+ ^qHy' (tf + iqHf it) f" it) - ^qf" it)' 



and 

Pdiff it) 



^(l-^-^«/W)/'W'-A/(t) (83) 
§,qH^f it)' + ^qHf it) f" it) - ^qf" it)' + ^qf it) J^^) (t) . 



Note that the equation of motion for the diff field depends on both H and H', 
where as the energy density and the pressure only depend on H. The equation 
of state, the pressure-to-energy density ratio, w it), for the diff field is given by 



wit) 



Pdiii it) 

Pdm it) 



(84) 



To study the properties of the diff field in the early universe, we assume 
that the contributions from the other sources arc negligible. In other words, we 
set the energy densities and pressures for other matters (including radiation) to 
zero. Under this assumption the Priedmann equations can be written in terms 
of the scalar field / {t) and its derivatives in time: 



H\t) = 



1 



/3 



A/(t) 



_ +^qH^f" it) + IqHf it) f" it) - ^qf" itf 



(85) 
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H'{t) 



(86) 



The expression of Hubble parameter H (t) can be obtained from Eq|S3 Since 
it is in a quadratic form in H, we have two possibihties: 



H±{t) 



9<z/' {ty 



iqf (t)f" {t)± 



(87) 



where 

Q = 



9qV' {tff'itf 



- + 9,/' itf) f 4 - 4) /' (i)^ + / it) (32A + 3«/' (t)^) \ 

^ " '^'^^n +g/"(t)^-2,/'(t)/(3)(t) y 

The Friedmann equations, Eq|SS| and Eg 1861 also allow us to simplify the 
equation of motion for Ea l81l these two equations can be used to eliminate 
H and H' in the equation of motion for Z?^^. 



5 Numerical Results 

We will numerically solve the equation of motion for the trace part of D^n in 
the Robertson- Walker metric. We use the numerical solution of / [t) to plot 
the Hubble parameter H (t), energy-density p (t), pressure p {t), the equation of 
state w (t) and p (t) + 3p (t) . Since the equation of motion that we are studying 
is very complicated, we will limit our analysis to the basic behavior of each 
quantities under given set of parameters and initial conditions. Our main goal 
in this numerical computation is not to obtain the results that agree with the 
current observation, but to investigate the effects of the trace part of diff field 
D^^ on the quantities mentioned above to study the nature of the scalar field 
/ (t). We use the explicit Runge-Kutta method for this numerical computation. 

5.1 Parameters 

EalSTI suggests that the Hubble parameter H (t) is expected to increase its 
magnitude very rapidly as 

fitf^-!^M^,. (89) 

The above condition also implies that the value of q should be negative. In this 
limit it is not difficult see that the positive solution (the solution with +\/Q) 
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Figure 2: The Trial Hubble Parameter Plot 



has a simpler behavior near the "peak" . Hence, in the following, we will study 
the positive solution case. Noting that the function of parameter P is basically 
to shift the value of kinetic term, we set this parameter to zero. Noting that the 
mass dimension scalar field / (t) is one, so it's natural to use the planck mass 
Mpi for this scalar field, i.e., we set Mpi to one. Hence undetermined parameters 
are a, q, and A. 



5.2 Trial Initial Conditions 

The equation of motion for the trace part of D^^i, in the Robertson- Walker metric 
is a fourth order nonlinear differential equation with three parameters (after 
setting (3 = and Mpi = 1). To solve these differential equations numerically, we 
need to specify trial values for the three parameters and four initial conditions. 
In this numerical computation the initial conditions for / (i), /' (i), /" (t), and 
/(■^^ (t) are used. A priori, we do not know how to choose these parameters 
and initial conditions. To overcome this, we use the physical condition on the 
Hubble parameter as a guide for choosing the trial values. We search for a set 
of the parameters and the initial conditions that give a positive and finite value 
for the Hubble parameter. With this as a criterion, a choice for the trial of 
parameters and the initial conditions become 

a = -l, q = ~6, A = 0.05, 
/ (0) = -10, /' (0) = 0, /" (0) = -1, = -1. (90) 

This set of trial values give the Hubble parameter plot shown in Fig[21 Note 
that we used k = (flat Universe), A = (no cosmological constant), /3 = 0, 
and A'Ipi = 1. 
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5.3 Computation 

As we have noted in Chapter 1, the nature of scalar fields are characterized by 
the equation of state w (t), and this value is defined as the ratio of the pressure 
p{t) to the energy density p{t). The acceleration of Universe can be studied 
through the quantity p{t) + 3p{t). We, thus, compute f (t), H{t), p{t), p{t), 
w {t) , and p {t) + 3p (t) for the reference set and study the effect of variations 
in two parameters {q and A) and four initial conditions of the scalar field / (t) 
if (0), /' (0), /" (0), and (0)) on these quantities. The following is the sets 
of two parameters {q and A) and four initial conditions of the scalar field / (t) 
(/ (0), /' (0), /" (0), and f^^^ (0)) used in this numerical computation 

5.4 Results 

The effect of changes in two parameters {q and A) and four initial conditions 
on the scalar field f (t) (/(O), /' (0), /"(O), and /f^) (0)) are studied. Each 
parameter/initial condition is varied with respect to the corresponding starting 
value. The results are plotted for the scalar field / (t), Hubble parameter H (t), 
energy density p (t), pressure P (t), and the equation of state (pressure-to-energy 
density ratio) w (t) . The unspecified values are the same as that of the reference 
set Eql^ni Here are our results. 

5.5 Effect of Changes in q 

Red (• • •) q = -2 , Green (— ) g = -6 , Blue (--)<? = -8 



f (t) 

-10 



-10 . 5 
-11 

-11.5 
-12 

-12 .5 



-13 . 5 
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Figure 3: Scalar field f{t) obtained with q = —2 (red), —6 (green), —8 (blue) 
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HCt) 



3 . 575 



0.55 



D.525 




Figure 4: Hubble parameter H{t) obtained with q = —2 (red), —6 (green) 
(blue) 




Figure 5: Energy density p{t) with q = —2 (red), —6 (green), —8 (blue 
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P[t] 




Figure 6: Pressure p{t) obtained with q = —2 (red), —6 (green), —8 (blue 




Figure 7: Equation of State w(i) obtained with q = —2 (red), —6 (green), 
(blue) 
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P<X) +3P(t) 




t 



Figure 8: p{t) + 3p{t) obtained with q = —2 (red), —6 (green), —8 (blue) 

5.6 Effect of Change in A 

Red (• • •) A 0.08 , Green (— ) A = 0.05 , Blue {- -) X ^ 0.02 

f (t) 




Figure 9: f{t) obtained with A 0.08 (red), 0.05 (green), 0.02 (blue) 



23 



HCt!) 

1. 65 

0. 6 
3 . 55 

0.5 
J . 45 

0.4 
2.351 




V. 



— T3" e— 10 



Figure 10: H{t) obtained with A = 0.08 (red), 0.05 (green), 0.02 (blue) 
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Figure 11: p{t) obtained with A = 0.08 (red), 0.05 (green), 0.02 (blue) 
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P[t] 




Figure 12: p{t) obtained with A = 0.08 (red), 0.05 (green), 0.02 (blue) 



^ 2 1 6 8 10 

\ 




Figure 13: Equation of State w{t) obtained with A = 0.08 (red), 0.05 (green), 
0.02 (blue) 
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P<X) +3P(t) 




Figure 14: p{t) + 3p{t) obtained with A = 0.08 (red), 0.05 (green), 0.02 (bl- 

5.7 Effect of Change in Initial Value of /(to) 
Red (• • •) f{to) - -8 , Green (— ) f{to) = -10 , Blue (- -) f{to) = -12 




Figure 15: Scalar field f{t) obtained with /(io) = —8 (red), —10 (green), — 
(blue) 
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Figure 16: Hubble parameter H{t) obtained with /(<o) = ^8 (red), —10 (green), 
-12 (blue) 



o.s 
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Figure 17: Energy Density p{t) obtained with /(to) 
-12 (blue) 



-8 (red), -10 (green), 
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P[t] 




Figure 18: Pressure p{t) obtained with /(io) = —8 (red), —10 (green), —12 
(blue) 




Figure 19: Equation of State w{t) obtained with f{to) — —8 (red), —10 (green), 
-12 (blue) 
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P<X) +3P(t) 




Figure 20: p{t) + 3p{t) obtained with f{to) = -8 (red), -10 (green), -12 (blue) 

5.8 Effect of Change in Initial Value of /' (to) 
Red (• • •) /' (to) = -0.5 , Green (— ) /' (to) = , Blue (- -) /' (to) = 0.5 

f (t) 




2 ■^4- — i-i-m ^_ 10 



Figure 21: Scalar field /(t) obtained with /'(to) = -0.52 (red), (green), 0.5 
(blue) 
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Figure 22: Scalar field H{t) obtained with /'(to) = -0.52 (red), (green), 0.5 
(blue) 




Figure 23: Energy density p{t) obtained with /'(to) = —0.52 (red), (green), 
0.5 (blue) 
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P[t] 




Figure 24: Pressure p{t) obtained with /'(io) = —0.52 (red), (green), 0.5 
(blue) 




Figure 25: Equation of state w{t) obtained with f'{to) — —0.52 (red), (green), 
0.5 (blue) 
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P<X) +3P(t) 




Figure 26: p{t) + 3p{t) obtained witli f{to) = -0.52 (red), (green), 0.5 (blue) 

5.9 Effect of Change in Initial Value of /" (to) 
Red (• • •) /"(to) = -0.5 , Green (— ) /"(to) = -1 , Blue (- -) /"(to) = -1.5 




Figure 27: Scalar field /(t) with /"(to) = -8 (red), -10 (green), -12 (blue) 



5.10 Effect of Change in Initial Value of (to) 
Red (• • •) /(3) = -0.5 , Green (— ) /(3) = -1 , Blue (- -) /(3) = -1.5 



32 



HCt!) 




Figure 28: Hubble parameter H{t) with /"(to) = -8 (red), -10 (green), -12 
(blue) 




t 



Figure 29: Energy density p{t) with /"(to) = -8 (red), -10 (green), -12 (blue) 
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P[t] 



t 




Figure 30: Pressure with /"(io) = -8 (red), -10 (green), -12 (blue) 




Figure 31: Equation of state w{t) with /"(to) = —8 (red), —10 (green), — 
(blue) 
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P<X) +3P(t) 




t 



Figure 32: p{t) + 3p{t) with /"(to) = -8 (red), -10 (green), -12 (blue) 



f (t) 




Figure 33: Scalar field f{t) obtained with f^^Hto) = -0.5 (red), -1 (green), 
-1.5 (blue) 
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Figure 34: Hubble parameter H(t) obtained with /'■^-'(io) = —0.5 (red), 
(green), —1.5 (blue) 
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Figure 35: Energy density p{t) obtained with /^"^^(io) = —0.5 (red), —1 (green), 
-1.5 (blue) 
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Figure 36: Pressure obtained with f^^^{to) = —0.5 (red), —1 (green), —1.5 
(blue) 




Figure 37: Equation of state 'w{t) obtained with f^^^to) — —0.5 (red), — 1 
(green), —1.5 (blue) 
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Figure 38: p{t) + 3p{t) obtained with f'^^^to) = -0.5 (red), -1 (green), -1.5 
(blue) 



6 Conclusion 

From the plots obtained in the last section we find that the equation of state 
w (t) can take, at least, values between —2 and 0. Most interestingly it always 
converges to —1 after some time for all the initial conditions we used. It means 
that the trace part of the diff field D/^^i, acts like a cosmological constant A 
in later times. This scalar field spends most of the time before w = — 1 as a 
phantom dark energy. This is a time-dependent phantom dark energy and its 
effects on Universe will be different from the time-independent cases jlUllllllTl] . 
Since the equation of state w can be in the very early Universe, it behaves 
as a dark matter when it happens. The fact that the value of w never gets 
close to 1/3 is very desirable because otherwise the radiation-like behavior of 
the scalar field at early times would influence already well explained primordial 
nucleosynthesis [7|. 

The effects on the scalar field on the Hubble parameter are also notable. 
It is easy to see that there are expected peaks [Eq|H7]. Figure^ Figure [TBI 
Figure 123 and Figure [2HI show shoulder- like regions where the value of Hubble 
parameter stay almost constant, i.e., the universe is expanding exponentially 
over this region: 

= 1 = constant ^ a (i) - (91) 

This scalar field, hence, acts like an infiaton over these regions. The constant 
Hubble parameter in later time can be interpreted as the current expansion of 
the Universe due to the cosmological constant like behavior of the scalar field. 
The last important observation is that the value of p (t) + 3p (t) is almost always 
negative. That is, this scalar contributes to the acceleration of Universe for most 
of time. 
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Though our rcsuhs are quahtativc, they arc showing very interesting cos- 
mological impUcations of the diflFeomorphism field 0^^: the trace part of the 
diffeomorphism field -D^,^ can produce the inflation-like behavior of the Universe 
in early time and the accelerating Universe in later time qualitatively. It is also 
observed that the time-dependent phantom dark energy plays a crucial role in 
this model. Our results, however, lack the realistic prediction: e-foldings N, 



during "inflation" predicted by this model is about A/' ~ 2, whereas any realistic 
model for inflation requires iV ^ 70. This suggest that more research in the 
study the set of parameters and the initial values necessary to make sensible 
predictions. Furthermore the understanding of the time-dependent phantom 
dark energy is essential in future research. 




(92) 
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